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(57) ABSTRACT 

A sensor system for detecting and estimating concentrations 
of various gas or liquid analytes. In an embodiment, the 
resistances of a set of sensors are measured to provide a set of 
responses over time where the resistances are indicative of 
gas or liquid sorption, depending upon the sensors. A concen- 
tration vector for the analytes is estimated by satisfying a 
criterion of goodness using the set of responses. Other 
embodiments are described and claimed. 

9 Claims, 6 Drawing Sheets 
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SYSTEM FOR DETECTING AND 
ESTIMATING CONCENTRATIONS OF GAS 
OR LIQUID ANALYTES 

PRIORITY CLAIM 5 

This application claims the benefit of U.S. Provisional 
Application No. 60/861,617, filed 29 Nov. 2006. 

GOVERNMENT INTEREST 10 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (35 USC 202) in which the 
Contractor has elected to retain title. 15 

FIELD 

The present invention relates to gas or liquid sensors, and 
systems for estimating concentrations of gases and liquids of 20 
interest. 

BACKGROUND 

Gas detecting sensors have many applications. One such 25 
application is to detect the presence of a dangerous gas, or to 
detect a gas whose presence indicates a dangerous situation. 

For example, NASA (National Aeronautics and Space Asso- 
ciation) has sponsored research in detecting low concentra- 
tions of sulfur dioxide (S0 2 ) in closed environments, such as 30 
for example aboard the international space station or the 
space shuttle. S0 2 could be a probable breakdown product 
from leaking lithium-thionyl chloride batteries. S0 2 is a col- 
orless gas or liquid under pressure with a pungent odor. Inha- 
lation or exposure could have adverse effects on human 35 
health. 

In addition to detecting various gases, or liquids, of inter- 
est, it may also be of utility to provide estimates of the con- 
centrations of various gases or liquids of interest. For 
example, various industrial processes may require that vari- 40 
ous gases or liquids be present within some range of concen- 
tration levels. In other applications, some dangerous gases 
may be tolerated if their concentrations fall below some 
threshold level. 

45 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a gas or liquid sensor system according to 
an embodiment. 

FIG. 2 illustrates a gas sensor according to an embodiment. 50 

FIG. 3 illustrates resistance values and response patterns 
for a gas or liquid sensor according to an embodiment. 

FIG. 4 illustrates a flow diagram for operation of a gas or 
liquid sensor system according to an embodiment. 

FIG. 5 illustrates a flow diagram for estimating gas or 55 
liquid concentrations gas according to an embodiment. 

FIG. 6 illustrates a prior art computer system for imple- 
menting various modules in the flow diagrams of FIGS. 4 and 
5 . 

60 

DESCRIPTION OF EMBODIMENTS 

In the description that follows, the scope of the term “some 
embodiments” is not to be so limited as to mean more than 
one embodiment, but rather, the scope may include one 65 
embodiment, more than one embodiment, or perhaps all 
embodiments. 


2 

FIG. 1 illustrates a sensor system according to an embodi- 
ment. Sensor chamber 101 includes a set of sensors for sens- 
ing various gases or liquids. This set of sensors, which may be 
termed a sensor array, is labeled as 103 in FIG. 1 . The set of 
sensors may or may not be a contiguous set of sensors, 
although in the particular embodiment of FIG. 1 the set of 
sensors is illustrated as contiguous. 

When the sensor system is operating, air is pumped from 
the surroundings into sensor chamber 101 . The air is directed 
either through activated charcoal filter 102, which is put in 
line to provide clean air for baseline data, or though a dummy 
filter of glass beads, dummy filter 104 , which is put in line to 
provide a pressure drop similar to that due to charcoal filter 
102 . Solenoid valve 106 is programmed to open the path to 
charcoal filter 1 02 and provide clean airflow for a pre-selected 
period of time at selected time intervals; otherwise, the air is 
directed through filter 104 . When air enters sensing chamber 
101 , whether through dummy filter 104 or activated charcoal 
filer 102 , the resistance of each sensor is measured. By pro- 
cessing the changes in resistances, embodiments provide for 
detection of various gases or liquids and estimation of their 
concentrations. 

Sensor chamber 101 may also include other sensors. In the 
particular embodiment of FIG. 1 , sensor chamber 101 
includes humidity sensor 105 to measure humidity. As dis- 
cussed later, measurements may be corrected for humidity. 
For some embodiments, humidity sensor 105 may be housed 
in a module separate from that containing sensor chamber 
101 . 

Microcontroller and data acquisition module 107 controls 
the various components in the sensor system, measures resis- 
tance, and records the acquired measurements. Microcontrol- 
ler and data acquisition module 107 may also include analog- 
to-digital conversion. Bus 108 provides an interconnect to 
other external components, such as a computer to process the 
resistance measurements. 

FIG. 2 illustrates a gas sensor according to an embodiment, 
where film 202 is formed on substrate 204 . For some embodi- 
ments, film 202 may include a polymer or co-polymer. Film 
202 adsorbs the gas, gases, liquid, or liquids of interest, which 
may be collectively referred to as the analyte. It is not known 
with certainty at this time if films according to the described 
embodiments absorb detectable analytes, adsorb detectable 
analytes, or some combination thereof. Accordingly, because 
the term sorption refers to the action of either absorption or 
adsorption, the term sorption and its variations will be used, 
indicating that a substance may be adsorbed, absorbed, or 
some combination thereof. 

Some examples of analytes are ammonia, mercury, sulfur 
dioxide, acetone, dichloromethane, ethanol, Freon 218, 
methanol, 2-propanol, toluene, and formaldehyde. Of course, 
this is merely a list of examples, and embodiments are not 
limited to detecting these analytes, and some embodiments 
may be designed without necessarily being able to detect 
these listed analytes. 

Substrate 204 may be flexible or rigid, and may be conduc- 
tive or non-conductive, depending upon the way in which the 
sorbed gas is detected. Sorption of the analyte of interest 
causes a change in one or more physical properties of film 
202 . By measuring this change, detection of a sorbed gas may 
be accomplished, provided the change is sufficiently large to 
allow a measurement. 

For example, for some embodiments, substrate 204 may be 
formed from silicon dioxide, and film 202 may include car- 
bon so that the electrical resistance of the carbon changes, 
depending upon the sorption of the analyte by film 202 . In this 
way, a resistance measuring device, in combination with the 
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system of FIG. 1, may provide for the detection of a sorbed 
analyte. Other embodiments may be based upon changes in 
other physical properties of film 202. For example, sorption 
of a gas increases the mass of the resulting system, so that a 
measuring device sensitive to mass change may be used in the 5 
detection of the sorbed gas. As another example, the sorption 
of gas by a co -polymer in film 202 may cause swelling, 
causing dispersion in an acoustic wave traveling through the 
system comprised of substrate 204 and film 202. Accordingly, 
an acoustic measuring system may be part of a gas sensor 
system. 

FIG. 3 illustrates the flow of signal processing applied to 
resistance measurements according to an embodiment of the 
present invention. For ease of discussion, the embodiment of 15 
FIG. 3 includes 32 sensors. The plots indicated by FIG. 3 need 
not actually be generated, but serve as representations of 
some of the signal processing functions performed by some 
embodiments. As indicated in plots 302, for each sensor a 
resistance is measured for some length of time. 20 

Plot 304 is a resistance vs. time plot for a particular sensor, 
say sensor number 32. Shown in plot 304 is a baseline resis- 
tance, denoted as R 0 . The baseline resistance is not necessar- 
ily measured directly, but may be estimated from resistance 
measurements. The resistance measurements may for some 25 
embodiments be obtained from raw resistance measurements 
after filtering has been applied. In the particular example of 
FIG. 3, an analyte has been introduced, and then removed, at 
various times, resulting in sudden changes in the measured 
resistance, R, as seen in plot 304. A time derivative of resis- 30 
tance, denoted as 

dR 

dt ’ 35 

is estimated from the resistance measurements, where t 
denoted the time index. (In practice, computations may be 
performed in finite arithmetic, so that mathematical opera- 
tions for some embodiments may be understood to be 40 
approximate with finite precision.) An event occurs when the 
time derivative 


dr 

is determined to be greater than some selected threshold, T. 
That is, an event is declared when 5Q 

dR/dt>T. 

For some embodiments, a baseline resistance associated 
with an event may be the measured resistance value at the 
beginning of the event. Accordingly, associated with each 55 
event is a baseline resistance. Baseline resistance may be 
denoted as R 0 . For some embodiments, the baseline resis- 
tance associated with an event may be taken at a time value 
other than the beginning of an event, or may be based upon 
filtered resistance measurements over some time interval 60 
associated with the event. 

When an event is declared, a differential resistance dR 
associated with an event is provided by taking a difference in 
resistances. For some embodiments, the differential resis- 
tance may be the difference between the resistance associated 65 
with an event taken at some time, say the beginning of the 
event, and the baseline resistance associated with that event. 
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For the particular example of FIG. 3, the baseline resis- 
tance is seen to be increasing over time, and the differential 
resistance dR is seen to be decreasing. Baseline drift may be 
due to various causes, such as changes in temperature, humid- 
ity, flow pressure, etc. For some embodiments, the baseline 
drift may not look like that of plot 304. 

The ratio of interest in plot 306 is dR/R 0 . The ratio dR/R 0 
may be referred to as a dynamic relative resistance change, 
where the term dynamic indicates that the baseline, as well as 
the differential in the resistance, may change with time. Plot 
306 illustrates the dynamic relative resistance change for the 
32 sensors, plotted as a bar graph. The ratio dR/R 0 may also be 
referred to as a response. 

In general, the set of responses dR/R 0 for all the sensors 
will yield a pattern indicative, to some degree, upon the vari- 
ous analytes present. This pattern may be referred to as a 
fingerprint. Based upon a fingerprint, embodiments estimate 
the presence and concentration of various analytes. For some 
embodiments, the responses dR/R 0 are adjusted for humidity 
by subtracting out the expected humidity response. For 
example, let dR(i)/R 0 (i) denote the response for sensor i, 
where i is an index ranging over sensor labels. Let dR^i)/ 
R 0// (i) denotes the expected, or estimated, response of sensor 
i due to humidity alone. Then, some embodiments may form 
a humidity -corrected response by subtracting out dR^ij/R^ 
(i) from dR(i)/R 0 (i) for each sensor i. 

The expected, or estimated, humidity response may be 
provided as follows for some embodiments. The set of 
responses may be measured in a controlled setting at various 
known humidity levels in which no detectable analytes are 
present so that a table of humidity responses for the various 
humidity levels may be generated and stored. In the field, the 
humidity level may be measured during real-time, such as for 
example with sensor 105 when air is being filtered by filter 
102 and the path through dummy filter 104 is closed. A 
look-up procedure applied to the stored table may yield an 
estimate of the expected set of humidity responses. 

As is well known, various interpolation and extrapolation 
techniques may be applied during a table look-up procedure. 
For example, in general the measured humidity level will not 
correspond exactly to any of the humidity level entries in the 
stored table, so that interpolation and extrapolation may be 
used. Other embodiments may correct the set of responses dR 
(i)/R 0 (i) for environmental factors other than humidity. For 
example, temperature may be taken into account. 

FIG. 4 illustrates a flow diagram according to an embodi- 
ment. Some or all of the modules indicated in FIG. 4 may be 
realized by dedicated hardware, software, firmware, or some 
combination thereof. Raw data is read from the output of the 
sensors (modules 402 and 404). Module 406 may include 
various signal processing functions, such as for example fil- 
tering to estimate baseline data, and processing to detect 
events. If an event is not detected (module 408), control is 
brought to module 404 so that additional raw data may be 
read. If an event is detected, then identification of target 
analytes is performed (module 410). 

If a known target analyte, or analytes, are detected (module 
414), the estimated concentrations are reported out and 
recorded (module 416). In the particular embodiment of FIG. 
4, if known target analytes are not detected, then various 
functional groups (families) may be identified (module 418). 
When testing for functional groups, subsets of the set of 
sensors are considered. That is, the fingerprint for a subset of 
the sensors is utilized in detecting whether a functional group 
is present. Different functional groups may correspond to 
different subsets of the sensors. Such functional groups may 
include, for example, an alcohol family (e.g., ethanol; 1-, 
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2-butanol; 1-, 2-propanol); an aromatic family (e.g., benzene; 
toluene; chlorobenzene; furan); a “one nitrogen” family (e.g., 
ammonia; acetonitrile; indole); and a ketone family (acetone; 
2-butanone). 

In the flow diagram of FIG. 4. note that control is brought 5 
to module 404 after module 408 if no event is detected, after 
module 418 if no functional group is detected, and after 
module 416 once a known target analyte or functional group 
is detected. However, such control flow is indicated merely 
for convenience, for in practice, various modules may be 10 
operating at the same time, or in a time-multiplexed fashion. 
Accordingly, the flow diagram indicated in FIG. 4 is not 
meant to imply any particular sequential order of operations. 

For example, module 404 may be operating in a continuous 
fashion, or in a near-continuous fashion if implemented in the 15 
digital domain, while other modules are operating. 

In determining whether a target analyte is present based 
upon a fingerprint of responses, for example where the fin- 
gerprint of responses may be the set of dynamic relative 
resistance changes 20 


where n is the number of sensors, some model of response 
versus the concentrations is chosen. This model may be rep- 
resented abstractly as 

7=jU*), 30 

where ~y is a vector of responses, x is a vector of concentra- 
tions, f (•,•) is a function, and A is a system model parameter 
describing system characteristics. In practice, A is a set of 
parameters, but for simplicity it may be referred to as a param- 35 

eter. The dimension of y is equal to or less than the number 
of sensors. For example, for some embodiments in which 
detection of a functional group is desired, the number of 
sensors under consideration may be a proper subset of the set ^ 

of sensors, so that the dimension of y is less than the number 
of sensors. When detection of one or more target analytes is 
desired, for some embodiments the responses of all the sen- 
sors are processed, so that the dimension of y is equal to the 
number of sensors. For some embodiments the dimension of 45 

x is equal to the number of target analytes. In practice, the 

dimension of x is less than or equal to the dimension of ~y. 

It has been found that for some embodiments, interesting 
results were obtained with the following model 50 

~y=A l x+A 2 x 2 , (1) 

where the superscript in x 2 is a notational convenience to 
mean that the \ th component of x 2 is the square of the \ th 55 

component of x , and A 1 and A 2 are matrices modeling the 
system characteristics. If the convention is chosen such that 

the dimension of y is n and the dimension of x is m, then for 
some embodiments using the above model, A x and A 2 are n by 6Q 
m matrices. 

The system modeling matrices A x and A 2 may be estimated 
by introducing known analytes of known concentrations, so 

that x and x 2 are known and y is measured. Standard 
numerical techniques may then be applied to Eq. ( 1 ) to find an 65 
estimate of the system modeling matrices. For example, as is 
well known, minimum norm least squares methods using 
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singular value decomposition may be applied to Eq. (1 ) to find 

estimates of A x and A 2 when x , x 2 , and y are given. 

Once the system modeling matrices are estimated, Eq. (1) 
then comes into play when an embodiment is deployed in the 

field, but where now A 1 andA 2 are known, y is measured, and 
x is to be estimated. Various well known numerical tech- 
niques may be applied to eq. (1) to estimate x for a given A : , 

A 2 , and "y , and one or more criteria of goodness are chosen 
for modules 414 and 418 to determine whether a target ana- 
lyte or functional group is present. 

It has been found for some embodiments that a numerical 
technique based upon a modification to the Levenberg Mar- 
quart Nonlinear Least Squares method, outlined in the flow 
diagram of FIG. 5, may be of utility in estimating the concen- 
trations from Eq. (1) for a given A 1? A 2 , and y . 

Referring to FIG. 5, the loop comprising modules 504 

through 516 updates a concentration vector x according to 
(module 516) 

X*-x+Ax , 

where Ax is a correction vector. 

A correction vector is generated as follows. Some initial 

starting point for the concentration vector 7, denoted as x 0 , 
is chosen in module 502. In module 504, each component in 

the residual vector ~y -f(A, 7) is weighted according to a 

weight vector w to provided a weighted residual vector r , 
where 

7=w*(y-J[A,'x), 

where * denotes component by component multiplication. 
The Jacobian matrix J is calculated in module 506, 

J=d r/dx, 

and the curvature matrix C is calculated in module 508, 

C=J f J, 

where f denotes transposition and complex conjugation (the 
Hermitian adjoint operation). In module 510, the matrix J w is 
calculated based upon the Jacobian and curvature matrices, 

J n =JC m , 

followed by a singular value decomposition 

J„=UDV t . 

A gradient vector g is calculated in module 512 according 
to 

In module 514, for some suitably chosen epsilon E, a correc- 
tion vector Axis calculated 

Ax = V(D 2 +eI)~ l/2 ~gC- l/2 , 

where I is an identity matrix chosen so that the dimension of 
€l matches that of D. 

The loop comprising modules 504 through 516 may be 
repeated for various values of €. For some embodiments, a 
series of candidate values for € such as €={ 10" 2 , 10 -1 , 1 , 10 2 , 
10 4 } may be chosen. On the first iteration of the loop com- 
prising modules 504 through 516, a search may be made 

through the series for € to find an x that results in the smallest 
norm (smallest sum of squares of the components) of the 
residual vector. For subsequent iterations, the previously used 
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values in the series for € may be multiplied by the previous 
found value for € that produced the minimum norm of the 

residual vector r . With the updated series, a search is again 
made for the best e as was done in the previous iteration. This 
process may be stopped after some chosen criterion of good- 
ness is satisfied. Criteria of goodness other than the norm of 
the residual vector may be considered in finding the best €. 

For example, a minimum norm criterion in conjunction with 
a criterion based on the least number of significant compo- 10 

nents in the x vector may be combined into a criterion of 
goodness. 

The above-described method for updating the concentra- 
tion vector over several iterations of the loop comprising 
modules 504 through 51 6 may then be performed for various 1 5 

starting points x 0 . For some embodiments, the various start- 
ing points may be chosen randomly (or more practically, 
pseudorandomly.) Some criterion of goodness, such as the 
norm of the residual vector, may be considered when choos- 20 
ing the final result for the various starting points. 

As discussed with respect to the flow diagram of FIG. 4, the 
various modules in FIG. 5 may be implemented in special 
purpose hardware, by software running on a programmable 
machine, firmware, or some combination thereof. FIG. 6 25 
illustrates, at a high level, a portion of a computer system, for 
implementing some of the modules illustrated in FIGS. 4 and 
5. For the computer system of FIG. 6, microprocessor 602 
communicates with memory 604 by way of chipset 606. 
Chipset 606 may be realized in one or more die, or some or all 30 
of its functions may be integrated with microprocessor 602. 
Memory 604 may represent a hierarchy of memory. Chipset 
606 also provides communication to I/O (Input/Output) mod- 
ule 608, which may be coupled to bus 108. Software for 35 
performing various modules in FIGS. 4 and 5 may reside in 
memory 604. 

Various polymers may be chosen for the sensor films. 
Examples are: Poly (4-vinylphenol-co-methyl methacrylate); 
Poly (ethylene-co -acrylic acid); Poly (styrene-co -maleic 40 
acid); Poly (2-vinyl pyridine); Poly (4-vinylpyridine); Vinyl 
alcohol/ vinyl butyral copolymer; Cyanoethyl hydroxy ethyl 
cellulose; Soluble polyimide; Polyepichlorohydrin; Poly (ep- 
ichlorohydrin-co-ethylene oxide); and Polyethylene oxide. 

Although the subject matter has been described in lan- 45 
guage specific to structural features and methodological acts, 
it is to be understood that the subject matter defined in the 
appended claims is not necessarily limited to the specific 
features or acts described above. Rather, the specific features 
and acts described above are disclosed as example forms of 50 
implementing the claims. Accordingly, various modifications 
may be made to the described embodiments without depart- 
ing from the scope of the invention as claimed below. 

Throughout the description of the embodiments, various 
mathematical relationships are used to describe relationships 55 
among one or more quantities. For example, a mathematical 
relationship or mathematical transformation may express a 
relationship by which a quantity is derived from one or more 
other quantities by way of various mathematical operations, 
such as addition, subtraction, multiplication, division, etc . Or, 60 
a mathematical relationship may indicate that a quantity is 
larger, smaller, or equal to another quantity. These relation- 
ships and transformations are in practice not satisfied exactly, 
and should therefore be interpreted as “designed for” rela- 
tionships and transformations. One of ordinary skill in the art 65 
may design various working embodiments to satisfy various 
mathematical relationships or transformations, but these rela- 
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tionships or transformations can only be met within the tol- 
erances of the technology available to the practitioner. 

Accordingly, in the following claims, it is to be understood 
that claimed mathematical relationships or transformations 
can in practice only be met within the tolerances or precision 
of the technology available to the practitioner, and that the 
scope of the claimed subject matter includes those embodi- 
ments that substantially satisfy the mathematical relation- 
ships or transformations so claimed. 

What is claimed is: 

1. A system comprising: 

a set of sensors, each sensor having a resistance indicative 
of gas or liquid sorption by the sensor; and 
a processing system to estimate values of a time derivative 
of resistance at various times for each sensor, 
the processing system to estimate event occurrences for 
each sensor by comparing for each sensor the estimated 
values of the time derivative of the resistance to a thresh- 
old, and the processing system to estimate baseline resis- 
tance for each sensor, and to provide for each sensor for 
each event occurrence a response indicative of a change 
in the resistance divided by the baseline resistance asso- 
ciated with each event occurrence. 

2. The system as set forth in claim 1, further comprising 
a sensor to measure humidity; 

the processing system to subtract from each response asso- 
ciated with each event occurrence a humidity response 
associated with each event occurrence based upon the 
measured humidity. 

3. A system comprising: 

a set of sensors, each sensor having a resistance indicative 
of gas or liquid sorption by the sensor; 
a processing system to estimate values of a time derivative 
of resistance at various times for each sensor; and 
a sensor to measure humidity; 

the processing system to subtract from a response associ- 
ated with an event occurrence for each sensor a humidity 
response associated with each event occurrence based 
upon the measured humidity. 

4. A system comprising, 

a set of n sensors, each sensor having a resistance R(i) 
indicative of gas or liquid sorption, where i E[1 , n] is an 
index indicating a sensor; and 
a processing system to process a set of resistance measure- 
ments for each sensor for a set of times, the processing 
system to estimate time derivatives dR(i)/dt for the set of 
times and to declare an event occurrence if there is a 
sensor i for which the estimate of the time derivative 
dR(i)/dt satisfies a relationship with a threshold; and 
the processing system to provide a set of responses dR(j)/ 
R 0 (j) for an event occurrence, where j is an index indi- 
cating a sensor belonging to the set of n sensors andj ES 
where S is a subset of [1 , n], S <= [1 , n] , where dR(j) is a 
difference in resistance values for sensor j associated 
with the event occurrence and R 0 G) is an estimated base- 
line resistance value for sensor j associated with the 
event occurrence. 

5. The system as set forth in claim 4, further comprising: 
a humidity sensor to provide a measure of humidity; and 
the processing system to provide a humidity response dR^ 

G)/Ro^) for j E S at the event occurrence based upon 
the measure of humidity, and to subtract the humidity 
response dR^j/Ro^) from the set of responses dRQ)/ 

R 0 G) to provide components for a vector response y . 

6. The system as set forth in claim 5, the processing system 

to provide a vector of concentrations x to satisfy a criterion of 
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goodness based upon a residual vector y (A, x ) where f(A, 

x ) is a function of the vector of concentrations x and system 
modeling parameter A. 

7 . The system as set forth in claim 6 , where the function 5 
f(A, x) is given by f(A, x)=A 1 x+A 2 x 2 , where the \ th 

component of x 2 is the square of the \ th component of x , 
where A l and A 2 are matrices representing system modeling 
parameter A. 

8. The system as set forth in claim 6, where the criterion of 10 
goodness is based upon a weighted residual 

7=W*(y_f(A, 7)), 

where * denotes component by component multiplication. 

9 . The system as set forth in claim 8 , where the vector of 

concentrations 7 is updated by an update vector Ax where 

Ax=V(D 2 +d)- 1/2 ~gC- 1/2 , 

where £ is a scalar, where a Jacobian matrix J is calculated, 


J=d 7/d7, 

and a curvature matrix C is calculated, 

c=7j, 

where f denotes transposition and complex conjugation, 
where a matrix J„ is calculated based upon the Jacobian and 
curvature matrices, 

J„=JC 1/2 , 

followed by a singular value decomposition 

J„=UDV t 

where g^is calculated according to 

‘g=Jn t 7. 

where I is an identity matrix chosen so that the dimension of 
El matches that of D. 



